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SUMMARY 

This report describes a preliminary investigation into the use of optimal control 
theory for the design of controls for a supersonic inlet. In particular, the control of a 
mixed-compression supersonic inlet is formulated as a linear stochastic optimal con- 
trol and estimation problem. An inlet can exhibit an undesirable instability (unstart) 
due to excessive inlet normal shock motion. For the optimal control formulation of the 
inlet program, a nonquadratic performance index is used, which is equal to the expected 
frequency of inlet unstarts. The computer program developed in this report determines 
optimal controllers which minimize this physically meaningful performance index. 

A linear lumped parameter, single-input, single-output inlet model is considered. 
The disturbance considered is a white Gaussian airflow perturbation at the compressor 
face station. The measured variable is normal shock position, corrupted by additive 
Gaussian white measurement noise. A Kalman filter is used to generate the estimates 
of the state variables. State estimate feedback gains are found by solving a series of 
quadratic control problems (Riccati equation) using the quadratic equivalence principle. 
The variation of the performance index with shock position tolerance, plant disturbance 
level, and measurement noise level is investigated. Results are presented for different 
levels of available control effort. 


INTRODUCTION 

For optimum performance of present day supersonic aircraft, it is necessary to 
provide an efficient means for decelerating the air from supersonic to subsonic velocity 
before it enters the compressor of the turbojet engine. A common method for accom- 
plishing this is a variable geometry supersonic inlet as shown in figure 1. The inlet 
shown is an axisymmetric mixed compression type to be used on future supersonic air- 
craft. In normal operation, air enters the inlet past a weak oblique shock wave and is 
compressed supersonically past the minimum area station (throat) up to the position of 
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Figure 1. - Schematic of supersonic axisymmetric inlet. 


the terminal normal shock. After the normal shock, the air is further compressed sub- 
sonically until it reaches the inlet of the turbojet engine (compressor face). For effi- 
cient operation, the inlet should maximize the recovered pressure at the compressor 
face over a wide range of operating conditions and maintain a stable flow pattern in the 
presence of external disturbances. 

The flow configuration of this type of inlet is bistable; the desired started condition 
just described and a second unstarted condition having a strong external shock system, 
subsonic flow in the initial convergent section, poor pressure recovery, and sometimes 
oscillatory flow (buzz). The desired started condition becomes unstable and switches to 
unstart wh^n either of two events occur: 

(1) A decrease in upstream Mach number causes the throat Mach number to decrease 
to one. In this case, a second (unstable) normal shock forms at the throat and moves 
rapidly upstream, out of the inlet. 

(2) A decrease in compressor face airflow causes the terminal normal shock to 
move upstream of the throat. Here again, the shock moves rapidly upstream, out of 
the inlet. 

The severe reduction in compressor face airflow and pressure occurring after an 
unstart may cause a compressor stall or combustor flame-out. In addition, the in- 
creased nacelle drag on an unstarted inlet may cause the aircraft to yaw suddenly. Also, 
a started inlet with normal shock too far downstream of the throat will have relatively 
poor pressure recovery and undesirable flow distortion at the engine face. Control is, 
therefore, required to maintain throat Mach number and shock position within acceptable 
limits to provide stable, efficient inlet operation. 

The primary inlet control variables are overboard bypass doors and a translating or 
collapsing inlet centerbody. Opening the bypass doors allows air to be dumped over- 
board, causing the shock to move downstream away from the throat. The translating 
centerbody varies the throat area, thereby varying the throat Mach number. A combin- 
ation of these two control modes is used to ensure stable inlet operation in the face of 
both upstream and downstream disturbances. 

The inlet can encounter both random and deterministic disturbances in a typical 
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flight environment. Deterministic disturbances are those such as atmospheric pressure 
changes due to a shock wave from a passing aircraft, pilot induced engine transients, 
aircraft maneuvers, and engine compressor stall. Disturbances which are random in 
nature are those such as atmospheric gusts, atmospheric temperature changes, combus- 
tion noise (in a duct burning fan) fed back from the engine, and aerodynamic compressor 
noise. 

Present inlet control systems (refs. 1 and 2) have been designed to minimize system 
response to deterministic disturbances such as ramp or sinusoidal changes in upstream 
Mach number and/or engine air flow rate. These systems were evaluated by examining 
the shape of their frequency response curves. Another approach which has been taken 
is to minimize system response to random disturbances. This was used by Barry 
(ref. 3) for the case where the random disturbance was atmospheric turbulence. He 
developed a method for predicting the expected number of unstarts per flight mile. Inlet 
controls were then evaluated on the basis of their effectiveness in minimizing unstart 
frequency. 

The approach described in this report is an extension of Barry's work to the prob- 
lem of designing optimal controllers which minimize unstart frequency. Basically, the 
inlet control problem is formulated as a linear stochastic optimal control problem using 
Barry's expected frequency of unstart equation as the performance index. This optimal 
design method is direct in that it leads to a control system which minimizes a physi- 
cally meaningful performance index rather than some intermediate quantity such as the 
area enclosed by the closed loop frequency response plot or shock position deviation. 

This performance index also leads to controllers tmiquely different from controllers 
which, for instance, minimize shock position deviation. 


OPTIMAL CONTROL FORMULATION 
Basic Formulation 

Bypass door control is used to avoid unstarts due to inlet disturbances which occur 
at the compressor face. Fast-acting bypass doors, such as described in references 
1 and 4, have been used to regulate against this type of disturbance. Also, to avoid un- 
starts due to atmospheric disturbances, the inlet throat area must be controlled. Typ- 
ically, this is done by translating the centerbody, which can only be done relatively 
slowly. However, the throat Mach number responds very rapidly to upstream (atmo- 
spheric) disturbances, due to the fact that flow is supersonic between the entrance to the 
inlet and the throat. Thus, it is difficult to mechanize a throat area control that will 
effectively reduce unstarts caused by upstream disturbances. For this reason, initial 
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Figure 2. - Single-input, single-output inlet control. 


studies have concentrated on the optimal control of bypass doors to reduce unstarts due 
to downstream disturbances only. 

Figure 2 is a block diagram of the inlet chosen for this study. Conventional control 
design for this inlet terminated with a long pipe ending in a choke plate was discussed 
in reference 1. For this report, an inlet transfer function model was obtained for the 
inlet terminated with a choke plate located at the compressor face station. It was found 
that, by curve fitting experimental frequency response data (unpublished NASA data), a 
good fit could be obtained over the range from 0 to 100 hertz using the form (Symbols 
are defined in the appendix. ) 



The dead time in G(s) is characteristic of inlet dynamics. However, in order to make 
use of the well-developed theory of linear stochastic optimal control and estimation 
(e.g. , ref. 5), a finite state (lumped parameter) model was needed. Thus, at the out- 
set, the dead time was approximated by a third-order Fade, resulting in a sixth-order 
state variable model for the inlet. 

In figure 2, the disturbance to the inlet w^^ represents a random air flow perturba- 
tion at the compressor face. It is assumed to be white Gaussian noise with a constant 
power spectral density (psd) equal to q. The output variable y_ (shock position) is 
measured through a noisy channel with measurement noise v . The noise is assumed 

o 

to be white Gaussian with a psd of r. 

The following performance index was chosen to be minimized; 

J = X + ko^ (2) 

% 
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where 



and 



expected frequency of inlet unstarts 
positive weighting factor 

mean square value of bypass door flow rate 
mean square value of shock velocity 
mean square value of shock position 


a shock position tolerance (distance between undisturbed shock position and inlet 
throat) 


J was selected so that the control must minimize unstarts X while limiting the amount 

of bypass door flow needed to do so. 

b 

The X relation (eq. (3)) is a classical exceedance equation given in reference 3. It 

gives the expected frequency with which the Gaussian random variable y_ exceeds the 

2 ° 

level a in the positive direction. Weighting factor k for ct is selected to penalize 

the control variable so that the level of control effort won't exceed that available. Selec- 
tion of the control effort weighting factor k will be discussed in a later section. 

In order to write equation (3) for the expected frequency of unstarts, the following 
assumptions must be made: (1) the inlet disturbances are Gaussian, (2) the inlet dynam- 
ics are linear, and (3) the controller is restricted to be linear and time- invariant. 

If a general approach to the optimal control problem were chosen (via dynamic pro- 
gramming, for example), the controller thereby obtained would, in general, be nonlinear 
due to the nonquadratic nature of the performance index of equation (2). This would be 
a violation of assumption (3) mentioned in the previous paragraph and, hence, the com- 
putation of X by equation (3) would be invalid. One way around this problem is to re- 
strict the controller to be linear and find the optimum linear, time -invariant controller 


which minimizes the expected frequency of unstarts. This is essentially what was done 
in this study. 
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It is known that the solution to the linear stochastic optimal control problem for 
quadratic performance indices dictates linear feedback of the optimal estimates of the 
state variables. This fact is a result of the so-called separation theorem (ref. 5). The 
solution technique can be extended to systems with nonquadratic performance indices by 
employing the quadratic equivalence principle developed by Edinger, Skelton, Stone, 
Ward, and White (ref. 6). The result is a linear control law which minimizes a quad- 
ratic approximation to the nonquadratic performance index. The method of quadratic 
equivalence was chosen for use in the inlet control problem because it forces a linear 
solution. In addition, the computer prc^rams required are those that have been fairly 
well developed for solving the linear optimal control problem for quadratic performance 
indices. 


Linear Stochastic Optimal Control and Estimation Problem Solution 

Before dealing with the nonquadratic problem, it is well to review the linear sto- 
chastic optimal control and estimation problem. It can be outlined as follows (ref. 5). 
Given a plant described by 


X = A(t)x + B(t)u + w (4) 

where x is the state vector, u is the control vector, and w is the plant disturbance 
vector. An output vector y is defined as 

y = H(t)x (5) 


and a measurement vector z is given as 

z = H(t)x + V (6) 


where v is the measurement noise vector. Both w and v are white zero mean Gaus- 
sian, and uncorrelated with each other, having correlation matrices defined by 


E|w(t)w'^(t + T I = Q(t)6(T) C7a) 

E|v(t)v'^(t + t)| = R(t)6(r) (7b) 
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A(t), B(t), and H(t) are real, possibly time varying, matrices. The initial state vector 
x(0) is assumed to be random with a covariance 

e{x(0)x'^(0;} = Pq (8) 

The problem is to minimize the quadratic performance index 

C = e|^- x’^(tj)S^^x(tj) + 1 j^x’^Q^(t)x + 2x"^N(t)u + u'^Pp(t)i:jdt| (9) 

Matrices S^, Qg(t), and N(t) are at least positive semidefinite, and Pg{t) is a positive 
definite matrix. 

The known solution to this problem (using the separation theorem) is 

u = -K^(t)x (10) 

where x is the optimal estimate of the state vector x and is generated with a Kalman 
filter described by 


X = A(t)x + B(t)u + Kg(t)[z - H(t)x] 


( 11 ) 


The optimal control gains K^(t) are given by 

Kc(t) = P’^t) [B’^(t)S(t) + n’^(I)] (12) 

where S(t) is the positive definite symmetric solution to the following matrix Riccati 
equation: 


-S - s(a - BP'^n"^) + (a - BPp - S^BP^^B’^js + 


S(tf) = S 



> (13) 


(For simplicity, the explicit time dependence of matrices A, B, etc. has been dropped. ) 
The Kalman filter gains K^(t) are given by 
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(14) 


Kg(t) = P(t)H(t)'^R(t)~l 


where P(t) is a positive definite symmetric solution to a Riccati equation, which is 


i> = AP + pa"^ - P(h’^R"^H)P + Q 

► 

P(0) = Pq 


(15) 


P(t) is also the covariance of the error in the estimate x. 

Equations (10) to (15) then define the solution to the problem of controlling a linear 
plant to minimize a quadratic performance index when the plant is subject to white noise 
disturbances. Figure 3 is a block diagram of the solution to the optimal control and es- 
timation problem showing the state estimator and state estimate feedback. The state es- 
timator (Kalman filter, eq. (11)) is basically a model of the plant driven by control u 
and measurement z. Signal z is compared with z, the estimated measurement, to 
form a term which is the error in the estimate of the measurement. This error is then 
multiplied by Kalman filter gains K^ and added back into the filter as a ’’correction" 
term. The, filter output x is weighted by the control gains K^(t) to form the optimal 


Plant noise w Measurement noise v 



Figure 3. - Combined optimal regulator-state estimator. 
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control vector u. In general, time varying control gain matrix K^(t) is calculated off- 
line by integrating equation (13) backward in time from t = tj to zero with the initial 
condition S(t£) = . Similarly, time varying filter gains Kg(t) are calculated off-line 

by integrating equation (15) forward in time, with the initial condition Pq being the 
a priori estimate of the covariance of x. 


Constant Gain Solution 

In selecting the nonquadratic performance index (eq. (2)), one assumes constant 
mean square values for all variables, this implying a constant gain system. That is, 
controller gains for the nonquadratic problem must be time invariant. For the quad- 
ratic problem, assuming matrices A, B, H, Q^, N, and P^ are constant, control 
gains Kp are constant for times sufficiently less than tj (which can be considered very 
large). is then given by 


=P"^(b'^S + n'^) (16) 

where S is the steady-state solution of the Riccati equation (eq. (13)). If, in addition, 

Q and R are constant matrices, Kalman filter gains become constant for time 
sufficiently greater than zero (so that initial condition effects have died out). is 
then given by 


Kg = Ph'^R"^ (17) 

where P is the steady-state solution of the Riccati equation (eq. (15)). The optimal 
system will then have all gains (Kg and K^) constant for times 0 « t « t^. During 
this time, all expected values are constant. 

Rewriting equation (9) as follows: 

C = e|- x"^(tj)SjX(t£)| + - J E^x'^QgX + 2x"^Nu + u'^Pgujdt (18a) 

it is seen that, for t^ — the major contribution to the index is the constant value of 
the integrand, obtained for t » 0 (after all initial condition effects have died out). 

Thus, the constant gain control derived can be said to have minimized that constant 
value given by 
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(18b) 


Cj = e|x"^Q^x + 2x"^Nu + u'^P^uj 

If the performance index for the inlet problem were purely quadratic, the solution 
to the control and estimation problem would be given by equations (16) and (17). How- 
ever, such is not the case. For the inlet control problem, a nonquadratic performance 
index defined by equation (2) is to be minimized. The following section will show how the 
principle of quadratic equivalence can be used to yield a solution which is a linear feed- 
back of state estimates identical to that of figure 3. 

Implementing Quadratic Equivalence 

For the quadratic problem, the necessary condition for the existence of a stationary 
point (such as a minimum) is that the variation of the performance index equal zero. 
First, define the scalars u 4 u^, w ^ W|^, and y ^ yg. Then, taking the variation of 
equation (18b) and setting it equal to zero, the following equation can be obtained: 

6e[x'^Q^x| + 26e{x%u^J + 6e[u^P^uJ = 0 (19) 

The variation of the nonquadratic performance index can be put in the form of equa- 
tion (19) by use of the quadratic equivalence principle. Then the known solution for the 
infinite time, linear, quadratic problem just discussed can be adopted. This can be done 
as follows. 

First, the variation of J (eq. (2)) is set equal to zero; that is. 



2 2 2 

Then equation (20) is expanded in terms of a , a. , and a giving 

^s ^s D 



where and W 2 are defined as 


( 21 ) 
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W« = 47rkcT a. exp/— 

^ ^ 8^8 I - ^ 



(22a) 


(22b) 


In order to compare equation (21) with equation (19), it is necessary to first ex- 
2 2 

press a- and a in terms of x, y, u, A, B, and H by using equations (4) and (5). 
^s ^s 

The mean square value of the output y written as an expected value, is 

s 




(23) 


but 


ys = Hx 


(24) 


thus. 


Oy = e|x’^h’^Hx| 


(25) 


Since 


Yo = Hx 


(26) 


then 


E|y2|= e[(Hx)’^(Hx)| 


(27) 


Thus, the mean square value of y becomes 

s 


a| = E{yg| = e{(HAx + HBu^ + Hw^)'^(HAx + HBu,^ + Hw^j)| (28) 
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The nature of the plant must be such that the product HW|^ is equal to zero; other- 
wise, CT. would be infinite and, hence, it would be impossible to minimize J. (This 
"s 

requirement is satisfied if the plant transfer function has more poles than zeros. ) With 


= e{, 


fTi rri rpi rri rp I 

k^AH^HAx + 2x^A^H^HBu^ + u^B^H^HBu^j 


(29) 


since Uv, is a scalar. 

b 2 

Also, a , as an expected value, can be written 


= ^K“b} = ^{“b%} <“> 

2 2 2 

Now, cr„ , , and a , from equations (25), (29), and (30) can be substituted into 

^s ^S % 

equation (21) to obtain 


f rp rp rp rp rp rp rp rp 

x^A^HAx + 2x^A^H^HBu^ + u^B^H^HBu^^j 


} + W^6e|x'^h'^HxJ 


+ W 26 E^ 




(31) 


Collecting terms results in equation (31) becoming 

6e{x'^(a'^h'^HA + WjH’^H^x} + 6e|2x'^^a'^h'^Hb)u^| + 6e|uJ(b'^h'^HB + W2)u^} = 0 

(32) 


Compare term by term equation (32), obtained from the nonquadratic performance 
index J, with equation (19), obtained from the quadratic performance C^. The following 
relations can then be made: 


= a'^h'^HA + 

N = a'^h'^HB 

p„ = b'^h'^hb + W, 

C Z 






(33) 
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The nonquadratic problem can now be solved by solving the quadratic problem (eqs. 
(13) and (15)) with Q^, N, and defined by equation (33), and P and S set equal to 
zero. Equivalence, however, only is assured at the optimum. At this point and 

Wq must satisfy equations (22a) and (22b) such that quadratic equivalence is maintained. 

^ 2 2 
Since and are defined in terms of mean square values such as cr and a- , 

^s 

the technique by which these values are determined is discussed in the next section. 


Mean Square Values 


2 2 

The mean square values a and a- are functions of P, the covariance of the 

"s ^s 

error in the estimate of the states, and X, the covariance of the states. Since the 
error in the estimate is 


e ^ X - X 


(34) 


the covariance of x can be written as 


X = e|xx'^| = E |(x - e)(x - e)'^J = X + P 
since plxe'^j = 0. Similarly, 


X = X + P 


(35) 


(36) 


where X is derived in Bryson (ref. 5, p. 417) as 


X = (A - BK^)X + X(A - BK^)"^ + 


(37) 


Using this value for X and obtaining P from equation (15), X becomes 


X = (BKj,)P + P(BK^)'^ + Q + (A - BKp)X + X(A - BK^)"^ (38) 

For the inlet problem, only the steady-state value of X is needed; thus, X = 0 and 
equation (38) becomes 


(BK^)P + P(BKj,)'^ + Q + (A - BK^)X + X(A - BK^,)'^ = 0 (39) 
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the steady- state covariance matrix equation. 

2 2 

Now, a and a- must be expressed in terms of P and X. From equation (25), 
^s 


Oy = e|x'^h'^Hx j = E |hxx'^H^ I = he Ixx”^ 


= hxh'^ 


(40) 


Ii 1 similar manner, equation (28) can be rearranged in terms of eIxx"^ ', e|xu^|, and 


ct| = HAE {xx"^! + 2HAE {xu^ 


b'^h’^ + HBE(u 




B^hT 


(41) 


But it can be shown that 


e{xuT) = -(X - p)kJ 


(42) 


and 




(43) 


Substituting these results into equation (41), becomes 

a? = HAXa'^h'^ - 2HA(X - P)kTb'^h'^ + HBK.(X - P)kJb'^h'^ (44) 

Yc ^ c c 


Thus Wj and W 2 can be calculated as functions of X and P using equations (40) 
and (44). In addition, mean square control effort is a simple function of X, given 
by' equation (43). 

This completes the derivation of the equations necessary to design an optimal inlet 
control system to minimize unstarts. Next, the solution of these equations will be de- 
scribed. 


CONTROL COMPUTATION PROCEDURE 

First, the estimation problem was solved, yielding a set of estimator gains K^. 
Then the quadratic control problem was solved for a number of trial (W^,W 2 ) pairs. 
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The solution to each quadratic problem K was used with the previously determined 

C Q e\ c\ c\ 

to compute the mean square values of the state vector. From this, , and 

were calculated. For a range of k’s, J, X, W^, and W 2 were then calculated. All 
these data were scanned off-line to determine optimum points. This was done by first 
finding points of minimum J for each value of k. As a check, an alternate method was 
used. This was to find points for constant k where the assumed was sufficiently 
close to the calculated and the assumed W 2 was close to the calculated W 2 . The 
accuracy with which the Wj and W 2 equations (eqs. (22a) and (22b)) were satisfied 



Figure 4. - Control computation procedure flow chart. 
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was a function of the number of trial pairs (Wj,W 2 ). That is, finding the optimum more 
accurately required more pairs. 

The digital computer program used in solving the inlet control design problem is 
outlined in figure 4. Inputs to the program are plant coefficient matrices A, B, and H, 
noise spectral densities Q and R, trial quadratic weights and W 2 , mean square 
control weight k, and output variable (shock position) tolerance a. 

The first step in the solution procedure is to compute the Kalman filter gains, which 
requires the steady-state solution of a Riccati equation. (Note that the Kalman filter 
gains need only be calculated once for a given pair of noise psd's. ) Next, the control 
weighting matrices Q^, N, and are found. Using these weights, the control prob- 
lem is then solved for the optimum feedback gains, where again a Riccati equation is 
solved. The mean square system behavior is then determined by solving the steady- 

state state covariance matrix equation. Then the mean square values of y, y, and u 

2 ,2 2 


can be found. The final calculations, using a. 


and , are W,, 
% ^ 


W 2 , A., and J, 


which are all printed out. 

The two major subroutines used are one for the steady- state Riccati equation and 
one for the steady- state state covariance matrix equation. The steady -state Riccati 
equation subroutine used the negative exponential method (ref. 7). For the inlet prob- 
lem, this method proved much faster than a direct fifth-order Runge-Kutta integration. 
The state covariance matrix equation was first solved by integration. However, since 
only the steady-state solution was needed, it was found that a faster method was to trans- 
form equation (39) to a set of n(n + l)/2 linear equations. Equation (39) is also a form 
of the so-called Lyapunov matrix equation. The solution method chosen was only one of 
a number of possible alternate methods for solving the Lyapunov equation (ref. 8). 


RESULTS AND DISCUSSION 
Introduction 

The preceding synthesis procedure was applied to develop an optimal estimator and 
control gains for the inlet described in equation (1). To do this, a sixth-order inlet 
model was obtained by combining a third-order dead time approximation with the third- 
order transfer function of equation (1). 

Parameters obtained from the curve fit using G(s) given by equation (1) are as fol- 
lows: 
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Kj = 5. 59 cm/(kg/sec) 
a = 210 rad/sec 


b = 80 rad/sec 
o)n = 365 rad/ sec 
C = 0. 30 
Tj = 0. 004 sec 

For a dead time of 0. 004 second, the third-order Fade is 

^-0.004 s _ -5.333x10'^° s^ + 1.60x10'^ s^ - 2.0xl0~^ s + 1 ^4, 

5. 333x10"^® s^ + 1.60x10’® s^ + 2.0x10’® s + 1 

Substituting for the dead time in equation (1) from equation (45) and using the aforemen- 
tioned constants, the sixth-order inlet transfer function becomes 


G(s) = 


Wj(s) 


(5. 59)(2. 001xl0^^)(-2. 537x10'^^ + 7. 083x10'^ - 7. 920xl0~^ + 2. 760xl0~^ s + 1) 

s® + 3. 30x10^ s® + 4. 80x10® + 3. 460x10® s® + 1. 160x10^® + 3. 230x10^^ s + 2. 001x10^® 


(46) 


Using equation (46), the A, B, and H matrices for the state variable representa- 
tion were conveniently obtained in phase variable form. But to improve the accuracy of 
the Riccati equation solutions, the problem was time scaled by 200 to 1, resulting in the 
following A, B, and H matrices: 


0.0 1.0 0.0 

0.0 — 1.0 


0.0 


A = 


1.0 


1.0 


1.0 


-312.7 -1009 -725.0 


-432.5 -120.0 


-16. 50j 
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= [ 0 , 0 , 0 , 0 , 0 , 1 . 0 ] 

H = { 5 . 59 [cm/(kg/sec)]| • [312.7, 172.6, -99.06, 17.71, -1.269, 0. 0] 
The plant disturbance vector was defined as 

= [ 0 , 0 , 0 , 0 , 0 , 


and its correlation matrix was 


E 


|w(t)w'^(t + t)| 


= Q6(t) 


where 


0 

0 

0 

0 

0 

0 


0 0 0 0 01 


qj 


and q is the psd of w^. Solutions to the estimation and control problem were obtained 
for a plant noise psd of 


q = 1.03x10"^ (kg/sec)^/(rad/sec) 
and two tvpical values of measurement noise v psd’s 

r - 3. 22x10"^ and 3. 22x10’® cm^(rad/sec) 

Note that, since w^ and v^ both have flat spectra (their psd’s not being functions 

of frequency), q and r do not change with time scaling. However, to compute X (un- 

2 2 

start frequency) in real time units, a and a • had to be expressed in real time. 

P ^s ^s 

Actually, only needed rescalir^. This is because 
”s 
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(47) 


(t) = E 


(ITlI 


and for r = yt, d/dt = y(d/dr), and (t) becomes 


= 4 ? ‘^ 3 1 


= r^E. 




2^2 

y ^,v 


(r) 


(48) 


2 2 
Thus, for y = 200, Oy (t) from the program was multiplied by (200) before X was 

calculated. ® 

A Kalman filter was designed for each q, r pair. Then a family of optimum con- 
trollers was designed by varying control weight k. The members of the family differ 
in that although all are optimum, each requires a different rms value of the control 

variable a, . Knowing a„ is important because it dictates how large the control 

% % 

actuator capacity must be. Conversely, if the maximum available control variable level 
is known, the controller can be designed to avoid control saturation. Saturation is un- 
desirable first because, if it occurs, the system becomes nonlinear and the whole opti- 
mum controller design procedure is in question. In addition, reliability is decreased 
if control saturation occurs. For example, in the case of the inlet, high deceleration 
forces are imposed on the bypass doors whenever they are driven full open or closed, 
which could lead to premature actuator failure. With these facts available, the results 
of the design program can be used to select a controller that will both be optimal and 
avoid saturation. 


Shock Position Tolerance Sensitivity 

The calculated optimal inlet control problem solutions can best be displayed by 
relating unstarts, control effort, and shock position tolerance for a family of optimal 
control systems. Figure 5 shows how well one family of optimum controllers is able 
to minimize unstarts. The expected frequency of unstarts X is plotted as a function of 
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rms bypass flow rate a for a typical range of shock position tolerance a. Each 

% 

value of a corresponds to a different controller (K„ gains), using the same Kalman 

filter (Kg gains) to generate the state estimates. Open loop values of X for different 
a's are shown along the vertical axis so that the improvement gained by optimal con- 
trol can be seen. Note that, as a goes to zero, the curves asymptotically approach 

% 

their open loop values, as would be expected. Also, as a increases, A decreases for 
both closed and open loop cases. 

A more conventional way of contrasting open and closed loop system performance 
is by comparing rms values. However, for the inlet problem, all disturbances have 
been assumed white and thus have infinite rms values. But if plant and measurement 
noise are assumed to be not white but colored, their mean square values become finite. 
Then various signals can be compared to them on an rms basis. For the following dis- 
cussion, assume that the spectral densities of W|^ (compressor face disturbance) and 
V- (shock position measurement noise) are flat to 2000 radians per second and fall off 

O 

on a Ic^-log plot at a two-to-one slope thereafter (first-order coloring). With such a 
high rolloff frequency, the noise will still appear essentially white to the inlet. The 
mean square value of W|^ can be calculated by integrating the psd of w^. For first- 
order coloring, mean square value of w^ turns out to be simply (ref. 9) 
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where q is the zero-frequency value of the psd and is the corner (rolloff) fre- 
quency of the coloring. 

The inlet being used has a nominal captured flow rate Wj of 16. 2 kilograms per 

second at Mach 2. 5. For a q of 1.03X10"^ (kg/sec)^/(rad/sec) and co^ of 2000 radians 

per second, the mean square value of w^^ is 1. 02 kilograms per second or 6. 3 percent 

of nominal inlet, captured flow. The open loop rms shock excursion for this disturbance 

level is then 1. 65 centimeters. For a low frequency measurement noise psd level of 

3.21x10" cm /(rad/sec) (shown in fig. 5), the rms noise-to- signal ratio for open loop 

shock position is about 11 percent. With control, rms shock position ranges from 

1. 65 centimeters down to 1. 14 centimeters for the maximum value of a shown. Al- 

% 

though this change is rather small, the corresponding change in > is large - about three 
orders of magnitude. This is because a (rms shock position) appears in the expo- 

nential term in A, thus A is quite sensitive to it. 

Note that as shock position tolerance increases, for a constant value of a , , the 

difference between closed and open loop A's increases. This indicates control is more 
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effective for large values of shock position tolerance. But larger A's cause lower inlet 

pressure recovery; thus there is a definite limit as to how large a can be. A trade-off 

then exists between A and a (pressure recovery). 

The other main control design consideration is what value of a is acceptable; 

% 

that is, how large can a be without having saturation occur. There is no completely 

% 

satisfactory answer to this question. The only sure solution is to simulate the system 

designed for a specific o , drive it with the desired random disturbances, and see 

% 

whether or not saturation occurs. Newton, Gould, and Kaiser (ref. 10) suggest, as a 

rule of thumb, o be no more than one-third of the saturation value to minimize the 

% 

possibility of control variable saturation. 

For the inlet studied, maximum bypass flow is about equal to nominal capture flow. 

Using Newton's rule, a suggested value of a would be 5. 4 kilc^rams per second. As 

% 

a design example, assume that it is desired that mean distance between unstarts be 

6 -7 

10 kilometers. For a cruise Mach number of 2. 5, this means that A must be 7. 5x10 

unstarts per second. Using figure 5, for a chosen to be 5.4 kilc^rams per second, 

"b 

shock position tolerance should be about 8 centimeters (see circular data point on fig. 5). 
Alternately, the inlet designer may insist some maximum value of o- not be exceeded 
so that sufficient pressure recovery is insured. In this case, figure 5 can be used to 
predict the expected distance between unstarts for some maximum and a. 


Measurement Noise Level Sensitivity 

In addition to shock position tolerance, other parameters which strongly influence 
the expected frequency of imstarts are plant noise psd q and measurement noise psd r. 
This follows from the fact that the ratio of q/r dictates the Kalman filter gains. The 
change in filter characteristics alters mean square values and thus A changes. Plant 
noise psd q is usually specified by the environment and cannot usually be changed by 
the designer. However, r is partly a function of the noise of the sensor and its asso- 
ciated instrumentation amplifier so the designer may be able to influence this variable. 

Figure 6 illustrates the effect a variation in r has on A. Plotted are three curves 
(solid lines) from figure 5 for cv = 6. 60, 7. 62, and 8. 64 centimeters plus curves (dashed 
lines) for lower measurement noise (where r is decreased by a factor of ten). The 
rms level of the low measurement noise is 3. 5 percent of open loop shock position rms 
compared to the 11 percent mentioned previously for the high measurement noise case. 
Measurement noise psd r has a strong effect on A, especially at higher rms control 
levels. For example, if, as in the previous example, a 10" -kilometer-per-unstart 
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figure is desired, it can be achieved with a of about 7. 4 centimeters instead of 8 cen- 
timeters if measurement noise is low. This means reliability can be held constant 
while pressure recovery is increased, if r can be reduced. Thus benefits can be 
gained by designing low noise transducers and amplifiers since they can lead to a 
smaller r. 


CONCLUDING REMARKS 

This report has described formulating the problem of controlling a supersonic inlet 
as an optimal control problem. A physically meaningful nonquadratic performance in- 
dex, the expected frequency of unstarts, was used in the formulation. A single-input, 
single-output inlet control case was solved via quadratic equivalence by solving a series 
of quadratic estimation-control problems. The method provides a straightforward pro- 
cedure for determining a set of linear feedback controls which minimize the nonquadra- 
tic performance index. Using this method for control design, the expected frequency of 
inlet unstart can be minimized. The effects of measurement noise and shock position 
tolerance were shown to have a large influence on the performance index. 

The theory used in this report can be extended to multivariable systems. Thus, 
assuming a suitable analytical model is available, control of both bypass doors and 
throat area, with downstream and upstream disturbances, can be included. Provisions 
can be made to allow coloring of both plant and measurement noise. The performance 
index can be expanded to allow for more than one mode of unstart. Also, more practical 
considerations can be introduced into the problem, such as how to minimize sensitivity 
of the state estimates to inaccuracies in the inlet model and prevent consequent degrada- 
tion in performance, inclusion of actuator dynamic limitations, use of other available 
measurements (such as compressor face pressure) to improve performance, techniques 
for simplifying the Kalman filter, and ways to handle dead time directly without intro- 
ducing Fade approximations. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 27, 1971, 

720-03. 
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APPENDIX - SYMBOLS 


A system matrix, n x n 

a transfer function zero, rad/sec 

B control matrix, n x c 

b transfer function pole, rad/sec 

C quadratic performance index 

time -invariant quadratic performance index 
c integer, number of control inputs 

e{ } expected value of 

e error in state vector estimate, n x 1 
G(s) inlet transfer function, cm/(kg/sec) 

H measurement matrix, m x n 

J nonquadratic performance index 

control gain matrix, c x n 
Kg estimator gain matrix, n x m 

Kj inlet gain, cm/(l^/sec) 

k positive weighting factor 

m integer, number of measurements 

N quadratic weighting matrix on state control, n x c 
n integer, number of state variables 

P solution of estimator matrix Riccati equation, n x n 
Pg quadratic weighting matrix on control, c x c 
Pq initial value of P, n x n 

Q power spectral density matrix of plant disturbance, n x n 

Qg quadratic weighting matrix on state, n x n 

q power spectral density of w^, (kg/sec) /(rad/sec) 

R power spectral density of measurement noise, m x m 

2 

r power spectral density of v , cm /(rad/sec) 

s 

S solution of control matrix Riccati equation, n x n 


25 


Sj terminal quadratic weighting matrix, n x n 

s Laplace variable, sec"^ 

t time 

tj final time 

u control vector, c x 1 

U|^ control bypass door flow rate, 1^/sec 

V measurement noise vector, m x 1 

Vg shock position measurement noise, cm 

Wj equivalence coefficient 

W 2 equivalence coefficient 

w plant disturbance vector, n x 1 

compressor face disturbance flow rate, kg/sec 
Wj total inlet flow rate, kg/sec 
X covariance matrix of x, n x n 

X cov 3 ,riance matrix of x, n x n 

X state vector, n x 1 

X estimate of state vector, n x 1 

y output vector, m x 1 

y_ shock position, cm 

z measurement vector, m x 1 

z estimated measurement vector, m x 1 

z_ measured shock position, cm 

a shock position tolerance, cm 

y time scaling factor 

6( ) first variation of 

6(t) delta function 

^ damping ratio 

X expected frequency of unstarts, unstarts/sec 

2 2 
a mean square bypass flow rate, (lig/sec) 

b 


26 


2 2 
a mean square disturbance flow rate, (kg/sec) 

2 2 
a mean square shock position, cm 

^s 

2 2 
a- mean square shock velocity, (cm/sec) 

^s 

T time, sec 

Tj inlet dead time, sec 

noise cutoff frequency, rad/sec 
natural frequency, rad/sec 
Superscripts: 

T matrix transpose 

differentiation with respect to time 
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